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1. INTRODUCTION 

Generation, transmission and distribution are the three stages of electrical power system. Generating 
power at generating station is transmitted through transmission lines and given to distribution network. 
It distributes power to loads connected to distribution network. Electrical distribution system is operated for 
lower voltages compared to transmission system. Because of lower voltage, higher currents pass in 
distribution system this lead to higher real power loss compared to generation and transmission systems. 
Distribution system is operated with fixed voltage and variable current magnitudes. In any environment three 
main important factors are existed. Those are security, reliability, efficiency (priority of choosing among 
three are security, reliability, efficiency). 

If any fault occurred in system, that may create cascade shutdown then overall system collapse. 
To avoid this should take immediate actions. Actions are may be isolated faulted zone in the system and give 
supply to remaining loads from healthy feeders, in some cases the remaining healthy feeders are already in 
overload condition then go for load shutdown at this moment system safety important role. This is nothing 
but service restoration, this may not provide 100% continuous power supply but somehow it maintains both 
reliable supply and system safety. 

The apparent power loss is equal to sum of real and reactive power losses. If power losses (apparent 
power loss) increases then real and reactive power losses also increases. Real power loss is more means 
system with less efficient. Reactive power is directly proportional to square of voltage in the system. 
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If reactive power is more, then voltage in the system more this lead to insulation breakdown on conductor 
finally dead short circuit will occur then overall system collapse, so need to maintain minimum apparent loss 
then it results good efficiency and voltage in the system. 

Power loss reduction formula is derived by Civanlar & Grainger [1], in this paper power loss 
reduction and reduced computation time done by using “branch interchange” algorithm. D. Shirmohammadi 
& H. W. Hong [2], in this paper the researcher focused on both resistive loss reduction and reduced 
computed time using DISTOP methodology. M.E. Baran & F.Wu [3], in this paper author applies two 
different load flows for the sake of accuracy and computation time. Both loss reduction and load balance 
achieved by using approximated power flow methods. Taylor, Lubkeman [4], developed a rule —based 
system aimed at the reduction of both search space and real power loss. J.Z. Zhu [5], presented to study 
distribution network reconfiguration based on refined genetic algorithm. A radiation distribution network 
load flow method is used to calculate power losses and objective is to minimize the system power loss. Jorge 
Mendoza et al [6], proposes a formal and robust strategy for applying the large optimal reconfiguration 
problems using genetic algorithm. The main objective is to power loss minimization through reconfiguration. 
The problems of overloading and service restoration, unbalanced load on electrical distribution system can be 
solved through reconfiguration by using ant colony optimization and heuristic search technique [7], [8] 
respectively. F.S Pereira, K. Vittori et al [9], proposed meta heuristic methodology for reconfiguration of 
electrical distribution system to minimize power loss. In this paper two algorithms ant colony optimization- 
travelling salesman and ant colony optimization-radial system reconfiguration are applied to one common 
standard five bus sample system. Finally author justified ant colony optimization-radial system 
reconfiguration is best for reconfiguration. Majid Jamil et al [10], proposed for the restoration of electrical 
distribution system in the overload condition with minimum power loss by using ant colony system 
algorithm. Firas M.F.Flaih,Lin, Lin Xiangning et al [11], presented the power loss minimization and voltage 
profile improvement in distribution system by using modified particle swarm optimization, modified particle 
swarm optimization was applied on the standard IEEE 33-bus test system. Other heuristic methods have been 
published that are principally based on switch selection(reconfiguration) done by using modified particle 
swarm optimization, Multi-Cooperative Particle Swarm Optimization, Selective Particle Swarm Optimization 
algorithms solved for IEEE 33 & 69 bus test systems. The objective function is to power loss minimization 
and approximate formulas for losses evaluation as those presented in [12], [13] and [14] respectively. 
S.Ganesh and R.Kanimozhi [15], proposed a new method for network reconfiguration in distribution system 
based on ant colony optimization with objective of real power loss minimization for 33-bus single feeder 
distribution network with five different kinds of load levels. Real loss minimization through reconfiguration 
of radial distribution system in both normal and abnormal conditions done with Improved Particle Swarm 
Optimization was projected by Dusharal Venkata Sunil et al [16]. 

In this paper, reconfiguration of distribution system has been formulated by using Artificial Ant 
Colony Optimization technique. Next output of any system with linear characteristics can be easily predicted 
but practically linear behaviour is not present. Nonlinear behaviour cases cannot predict next output, as it is 
very complex. This technique gives the solution to nonlinear problems, reduces time consumption and gives 
accurate results. This technique is co-related with electrical distribution system as current matches with ant’s 
movement, pheromone intensity matches with power flow and length of the line with resistance or time 
consumption. Two objectives are met, first one is service restoration and second one is real power loss 
minimization. The proposed work is organized into six sections, first Section-1 is introduction, Section 2 on 
problem formulation, Section 3 discusses proposed technique algorithm for reconfiguration, and Section 4 
discussed proposed method applications on test systems. Section 5 on results and finally conclusions are in 
Section 6. 


2. PROBLEM FORMULATION 
2.1. Service Restoration 

First identify faulted feeder, isolate the fault. Identify the affected zones due to fault, which are out 
of service. Reconfigure the system with LBI calculations using (1). 


LBI=-J¥W -¥,) (1) 


Where: n= no. of primary feeders; Y=average of the normalized loadings Y; 


i‘ —Feeder Actual Loadin i 
Y= rt ae Oe 


ith_Feeder Base Loading 
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Subjected to following constraints: 
Si < Sri 
Pi< Pri 
Qi Qari 
Where: i=1, 2, 3....n; n= no. of primary feeders 
S=MVA loading on feeder — i ; Spj=Rated MVA on feeder - i 
P;=MW loading on feeder — i; PR=Rated MW loading on feeder - i 
Q=MVAR loading on feeder —i ; Qg=Rated MVAR on feeder — i 
In each and every reconfiguration load should be in balance condition, radiality and critical 
parameters with in permissible limits. Among all possible combinations which one gives least LBI value and 
satisfies all critical parameters that is the optimal solution. 


2.2. Power Flow Analysis - Direct Load Flow Method (DLFM) 
This method performs the power flow analysis for weakly meshed distribution system under balanced 
operating condition employing constant load model. In this method three steps are present, those are 
a. Equivalent current injection 
b. Formation of BIBC matrix 
c. Formation of BCBV matrix 


2.2.1 Equivalent Current Injection 
It is purely based on current injection. At i bus, the apparent power S; is specified and the 
corresponding equivalent current injection at the n" iteration of the solution is computed as: 


S| =P; +jQi 2 
. Si_)« 
P= (eve) +5) = (Ga) 3) 


Where: i=1, 2, 3,4... N 

S; is the apparent power at i bus; P; is the real power at i" bus. 

Q; is the reactive power at i" bus; V? is the bus voltage at the n" iteration for i" bus. 

I? is the equivalent current injection at the n" iteration for i” bus. 

[Real & ]”are the real and imaginary parts of the equivalent current injection at the n" iteration for i” bus. 


2.2.2 Formation of BIBC Matrix 
In Figure | is used as an example to understand the formation of BIBC matrix. Using (3) apply 
Kirchhoff current law (KCL) to distribution network. 





Figure 1. A simple weakly meshed distribution system 


B1=12+13+14 (4) 
B2=13+14 (5) 
B3=14 (6) 


Above (4), (5), (6) can be written in matrix form. This matrix is nothing but BIBC (Bus Injection to 
Branch Current) matrix as shown below (7). 
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B, 11 1)fl 
B -| 1 i (7) 
B; 00 ULL 
[B] = [BIBC] uxcv-1) [1] (8) 


Where: 

Bl, B2, B3...... BN are the branch currents. 

V1, V2, V3.....VN are the bus voltages. 

Z12, Z23, Z34.....Z3N are the branch impedances between land2, 2 and 3, 3 and 4, 3 and N respectively. 
12, 13, I4......... IN are the load currents. 

M is the number of branches. 

N is the number of buses. 


2.2.3 Formation of BCBV Matrix 
BCBV matrix nothing but, it is the relation between branch currents and bus voltage. 


V1=V2 + (B1*Z12) (9) 
V2=V3 + (B2*Z23) (10) 
V3=V4 + (B3*Z34) (11) 


By using above (9), (10) and (11) the voltage at Bus-4 can be rewritten as 
V, = Vi — (ByZ12) — (B2Z23 ) — (B3Z34) (12) 


Above (9), (10), (11) and (12) can be written in matrix form. This matrix is nothing but BCBV 
(Branch current to Bus Voltage) matrix as shown below (13). 

















Vi V, Z42 0 O 7 /B, 

Vi} —|V3} = |Zi2 Zo3 0 le (13) 

Vy V4 Z42_ Z23 Z3411B3 

[AV] = [BCBV] (N-1)XM [B] (14) 
Pi+Qi |" 

p= PES (15) 


Update voltages using (14), Calculate DLF (direct load flow) matrix using (16) 
[AV]"** = [DLF][I]" (16) 


Set iteration count n=0, if |[I?** — I?]| < & (tolerance value) then calculate line flows and losses using final 
voltages. Otherwise go for iteration number increment (n=n+1) and follow above same procedure. 

Subject to the following constraints: 

a. Radiality to be maintained. 

b. Number of tie switches before and after reconfiguration should not change. 

c. Voltage Magnitude: 


Vmin </lV;l < Vinax 31 (17) 
Where: 


V; is the voltage on i™ bus; Ne is the number of branches in the system; N, is the number of buses in the 
system. These constraints are applicable to all Objectives. 
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3. PROPOSED TECHNIQUE - ARTIFICIAL ANT COLONY OPTIMIZATION (AACO) 
3.1. Introduction 

This technique gives solution to nonlinear problems related to any field of work. Natural living 
behavior of ants gives optimal solution either a minimal or a maximal of main objective of problem is called 
Ant Colony Optimization. 

When an ant moves to search for food from home to food source it deposits chemical substance on 
ground called as pheromone. This substance is a clue for other ants from the same home to collect the food. 
Ants follow the pheromone trails left by the first ant. From home to food source the ants have several paths 
with variable pheromone trails. Among those trails ants choose higher pheromone intensity path. This path is 
also attracted by remaining ants. Always ants prefer shortest path (i.e. higher pheromone intensity). 

Figures 2 and 3 represent ant colony system without and with obstacle. Thicker line indicates higher 
pheromone intensity as a shorter path, thinner line indicates less pheromone intensity as a longer path. 





Figure 2. Natural behavior of ant without obstacle Figure 3. Natural behavior of ant with obstacle 


Electrical Distribution System having n= number of switches and total possible combinations for 
reconfiguration are 2".In order to reduce time consumption for this possible reconfigurations ACO itself has 
probabilistic transition rule (PTR). 

PTR measures ant probability of choosing path among all possible paths. For example when a system 
has 10 switches then total number of possible combinations for reconfiguration are 2'°=1024. Using PTR this 
may reduce to 15 possible combinations for reconfiguration and gives optimal solution. 

PTR is formulated as (18) 


[tj (t)] “In ij (t)] R 


ust ) Eales)“ [nisl 


(18) 
Where: j, q& NK. 


NK is a list in the memory of ant that records the cities which will be visited to avoid stagnations. 
i, j is the travel path of ants; t;,;=pheromone trail deposited between the nodes i, j. 


n;,;==represents the heuristic information=_— ; dj j= distance between city i and j. 
ii” 


«, 8 =Are two parameters that influence the relative weight of pheromone trail and heuristic guide function 
respectively. 

After each tour is completed a local pheromone update is determined as below mentioned (19). It is 
the sum of old pheromone intensity (p. t;;(t)) and new pheromone intensity(At;;). 


Ty(t+ 1) = 4, + (At;;) (19) 


Where: 
-yl k 
Ati = pa Ati, 
Atk = 2. if ant k travels on edge (i, }) otherwise it is zero; t =iteration counter; k= 1, 2....1 ants. 
; k 


Q=constant; L,=length of tour. 
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3.2. Flow Chart 
Flowchart of proposed method is shown in Figure 4. 
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Figure 4. Flowchart of proposed method 





4. PROPOSED METHOD 
The proposed AACO method is tested on four feeder service restoration, Ohmic loss minimization 
of IEEE 33 bus and 69 bus test systems. Test results as follows below. 


4.1. IEEE Four Feeder Service Restoration 

The proposed AACO technique is tested on the IEEE four feeder test system is shown in Figure 5, 
in this 16 sectionalized switches are generally in closed condition(between zone loads example-switches 
2,3,4...) and 15 tie switches are generally open condition (between zone loads of different feeders example- 
switches 6,7,8...). Each feeder loading capability is 10,000 kVA.1, 11, 21, 31 are the circuit breakers of four 
feeders in the system. Table 1 shows the zone loads. 





Figure 5. Four feeder test system 
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Table 1. Four Feeder Load Data 








Zone number Zone load(k VA) 
1 1600 
2 700 
3 1800 
4 500 
5 1900 
6 500 
7 1500 
8 1000 
9 500 

10 800 
11 3000 
12 3500 
13 700 
14 1000 
15 600 
16 1500 
17 2000 
18 1800 
19 1500 
20 500 





From Figure 5, each feeder capacity is 10000 kVA and actual load on feeder A, B, C and D is 6500 
kVA, 4300 kVA, 8800 kVA, 7300 kVA respectively. For illustration of the proposed technique, assume 
some random faults occur on four feeder system (i.e. feeder A, B, C and D at positions are Z1, Z8, Z14 and 
Z17 respectively).these random faults explained in below four cases. 

Case-1: Fault on feeder A at switch-1 (circuit breaker) 

Because of fault at switch-1, faulted section is isolated by opening Switch-1(Circuit breaker).To 
maintain continuous power supply need to connect total (Or) at least some of critical loads are transfer to 
other healthy feeders. 

Before Reconfiguration: load on feeder A is zero kVA, feeder B is 4300 kVA, feeder C is 8800 kVA and 
feeder D is 7300 kVA. 

After Reconfiguration: Switches 6, 20 are closed condition (sectionalized Switches), Switch 15 is in open 
condition (tie switch).Load on feeder A is zero kVA, feeder B is10000 kVA, feeder C is 9600 kVA and 
feeder D is 7300 kVA. 

Case-2: Fault on feeder B at position Z8 

Because of fault, faulted section is separated by opening switches 13 and14. Discarded load on 
feeder B is 2300 kVA. 

Before Reconfiguration: load on feeder A is 6500 kVA, feeder B is 2000 kVA, feeder C is 8800 kVA and 
feeder D is 7300 kVA. 

After Reconfiguration: switch 9 is in closed condition (sectionalized switch). 

Load on feeder A is 7800 kVA, feeder B is 2000 kVA, feeder C is 8800 kVA and feeder D is 7300 kVA. 
Case-3: Fault on feeder C is at position Z13 

Because of fault, faulted section is separated by opening switches 23 and 24. Discarded load on 
feeder C is 2300 kVA. 

Before Reconfiguration: load on feeder A is 6500 kVA, feeder B is 4300 kVA, feeder C is 6500 kVA and 
feeder D is 7300 kVA. 

After Reconfiguration: switch 19,30 are in closed condition (tie switches),switch 35 open condition(i.e. 
Sectionalizing switch). Load on feeder A is 6500 kVA, feeder B is 6400 kVA, feeder C is 6500 kVA and 
feeder D is 6800 kVA. 

Case-4: Fault on feeder C is at position Z14 

Because of fault, faulted section is separated by opening switches 24 and 25. Discarded load on 
feeder C is 1600 kVA. 

Before Reconfiguration: load on feeder A is 6500 kVA, feeder B is 4300 kVA, feeder C is7200 kVA and 
feeder D is 7300 kVA. 

After Reconfiguration: switch 20 is in closed condition (sectionalized switch).Load on feeder A is 6500 kVA, 
feeder B is 4900 kVA, feeder C is 7200 kVA and feeder D is7300 kVA. 

Case-5: Fault on feeder D at position Z17 

Because of fault, faulted section is separated by opening switches 32 and 33. Discarded load on 
feeder B is 5800 kVA. 
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Before Reconfiguration: load on feeder A is 6500 kVA, feeder B is 4300 kVA, feeder C is 8800 kVA and 
feeder D is 1500 kVA. 

After Reconfiguration: switches 10, 19, 29 are in closed condition(Sectionalized switches) and switches 15, 
24 are in open condition (tie switches).Load on feeder A is 7800 kVA, feeder B is 2000 kVA, feeder C is 
8800 kVA and feeder D is 7300 kVA. 


4.2. Minimization of Ohmic Losses 
4.2.1 TEEE 33 Bus Test System 

The proposed technique was applied on IEEE 33 bus test system. The test system base voltage is 
12.66kV,base MVA -100 and total system real and reactive loads are 3715kW and 2300kVAR, respectively, 
and it consists of five tie switches. 

Before reconfiguration - Initial network configuration consists of S1 to S32 are the sectionalizing 
switches and S33 to S37 are the tie switches. For this combination real power, reactive power and apparent 
power losses are 202.677 1kW, 135.1410 kVAR and 243.6 kVA respectively. 

After reconfiguration - the optimal configuration obtained by the proposed algorithm is S29, $33, 
S34, S35, S37 are the tie switches (i.e. dotted lines) remaining all are the sectionalizing switches as shown in 
below Figure 6. For this combination real power, reactive power and apparent power losses are 109.6065kW, 
78.0381 KVAR and 134.5494 kVA respectively. After (Optimal) reconfiguration, 45.9206 % of loss reduction 
is obtained comparatively. 


1 
sl 
s18 2 
aa =e $22 
s19 s3 23 
20 i $23 
520 s4 Tai 
21 $33 5 <— 95 
2 : - : 
22 1 6 25 
1 1 oh Z 26 : 
I am ( s26_ 4 1837 
| 8 s27 1 
I s8 28 I 
I = oe $28 Do pes 
I 10 I 1829 
“aed s10 . I = 30 
L sll 1 —— 31 
= ee ee 2 1s34 s31 
$12 " af 
s32 
4 13 1 33 
s 
14 I 
sl —_ = of 
15 
15 
sey 
s16 536 


Figure 6. 33 Bus test system after reconfiguration 


4.2.2 TEEE 69 Bus Test System 

The proposed technique was applied on IEEE 69 bus test system. The test system base voltage is 
12.66 kVA, base MVA-100 and total system real and reactive loads are 3802.19kW and 2694.6kVAR, 
respectively, and it consists of five tie switches. 

Before reconfiguration-Initial network configuration consists of switches 1 to switch 68 are 
considered by the sectionalizing switches and switch 70 to switch 74 are the tie switches. For this combination 
real power, reactive power and apparent power losses are 225.0028kW, 102.1659 kKVAR and 247.1116 kVA 
respectively. 

After reconfiguration-switches 49, 64, 70, 71, 72 are the tie switches (i.e. dotted lines) remaining are 
the sectionalizing switches as shown in below Figure 7. For this combination real power, reactive power and 
apparent power losses are 71.1381kW, 78.7572 kVAR and 106.1288 kVA respectively. After (Optimal) 
reconfiguration, 68.3834 % of loss reduction is obtained comparatively. 
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Figure 7. 69 bus test system after reconfiguration 
5. RESULTS 


5.1. IEEE Four Feeder Service Restoration 
The simulation results of distribution system reconfiguration based on AACO for IEEE 4- feeder 
service restoration are shown in Table 2. 


Table 2. Results of 4 Feeder System under Faulty Conditions at Different Locations 








Feeder Fault at Optimal switching operation LBI 

A Switch-1 Closed Switches - 6&20 0.200853 
Open Switch - 15 

B Z8 Closed Switch - 9 0.131974 

Cc Z13 Closed Switches - 19&30 0.0075 
Open Switch - 35 

Cc Z14 Closed Switch - 20 0.048007 

D Z17 Closed Switches - 10,19&29 0.140506 


Open Switches - 15&24 





In IEEE four feeder test system, for same fault on feeder-C at position Z13. Continuity of power 
supply maintained, optimal LBI value of proposed method is compared with other reference papers LBI 
values, among those proposed method gives least LBI value as shown in below Table 3. 


Table 3. Results Validation of Four Feeder System for Service Restoration 








Method LBI value for the same fault on feeder-C at position Z13 
Dr.P.Ravi Babu[22] 0.0100 
Dr.P.Ravi Babu[23] 0.0124 
Proposed method 0.0075 





5.2. Minimization of Ohmic Losses 

In order to express the proposed method IEEE 33and 69 bus test systems are taken into 
consideration. The method is coded using MATLAB version R2017a is used to run the optimal power flow 
solver with direct load flow method. AACO is carried out best results, are shown in Table 4 and Table 5. 


Table 4. Results Validation for 33 Bus Test System 








Method Tie Switches Loss After Reconfiguration(kW) %I-R Loss Reduction 
J.Z.Zhu [5] $7,S9,S14,S32,S33 139.532 31.1555 
Majid Jamil et al. [10] $2,$13,816,S33,S37 112.68 44.4041 
Firas M.F.Flaih et al. [11] $7,89,S14, $32,S37 139.55 31.1466 
proposed $29,S33,S34,835,S37 109.6065 45.9206 
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Table 5. Results Validation for 69 Bus Test System 








Method Tie Switches Loss After Reconfiguration(kW) %I°R Loss Reduction 
Inji Ibrahim atteya et al. [12] 14,55,61,69,70 100.6 55.2894 
Soumitri Jena et al. [13] 12,18,58,61,69 103.62 53.9472 
Tandon.A et al.[14] 13,20,55,61,69 107.05 52.4228 
proposed 49,64,70,71,72 71.1381 68.3834 





Improvement in Bus Voltages Before and After Reconfiguration of 33 Bus and 69 Bus Test 


Systems, as shown in Figure 8. 


Before and After Reconfiguration Bus Voltages in Before and After Reconfiguration Bus Voltages in 
33 Bus Test System 69 Bus Test System 


123 4 5 6 7 8 9 101112131415161718 1920212223 24252627282930 313233 
1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 61 64 67 


—— Before Reconfiguration Bus Voltage(Kv) —— After Reconfiguration Bus Voltage(Kv) —— Before Reconfiguration Voltage(kv) ——After Reconfiguration Voltage (Kv) 





6. 


Figure 8. Bus voltages before and after reconfiguration of 33 bus and 69 bus test systems 


CONCLUSION 
Service restoration and Ohmic loss minimization in distribution system are the main objectives in 


this paper. Service restoration achieved on IEEE 4-feeder test system under faulty conditions at different 
locations and for the same fault at Z13 on feeder-C least LBI value is obtained with proposed technique 
compared to other reference papers. The Ohmic losses are minimized optimally in IEEE 33 and IEEE 69 bus 
test systems with I? R loss reduction by 45.9206% and 68.3834% respectively. The minimum voltage before 
reconfiguration is 11.5597 kV at node-18 and after reconfiguration it is raised to 11.7876 kV at node-18 for 
33 bus test system. The minimum voltage before reconfiguration is 11.5103 kV at node-54 and after 
reconfiguration it is raised to 12.181 kV at node-69 for 69 bus test system. These objectives achieved through 
reconfiguration by using artificial ant colony optimization (AACO). 
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